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Abstract. The radio light curve and spectral evolution of the blazar CTA 102 during its 2006 outburst can be 
rather well explained by the standard shock-in-jet model. The results of a pixel-to-pixel spectral analysis of multi- 
frequency VLBI images, together with kinematics derived from the MOJAVE survey lead to the picture of an 
over-pressured jet with respect to the ambient medium. The interaction of a traveling shock wave with a standing 
one is a possible scenario which could explain the observed spectral behavior. 



1. Introduction 

The bl azar CTA 102, (B2230+11 4), has a redshift of z = 
1.037 (|Hewitt fc Burbided Il989h and it is classified as 
a high polarization qu asar (HPQ) with a linear o ptical 
polarization above 3% (jVeron-Cettv &: Veronl 120031 ) with 
an optical magnitude of 17.33. The source wa s observed 
for the first time bv lHarris fc Roberts! ( 1960h and right 
from the beginnin g it was showing radio variability (e.g., 
Sholomitskii 19651 ) which led other scientists to suggest 
that the signal was comin g from an extraterrestrial civ- 
ilization (jKardashevI Il964f) . Later observations identified 
CTA 102 as a quasar. 

Since that time CTA 102 has been the target for nu- 
merous observations at different wavelengths. Besides the 
mentioned variation in the radio flux density, CTA 102 
changes its optic al behavior in a no less spectacular way. 
Pica et al. ( 1988f) reported rapid variation up to 1.14 mag 
around an average value of 17.66 mag. An increase of 
1.04 mag within 2 days in 1978 was so far the most signif- 
icant outburst. The strongest radio flare since 1986 and 
a nearly simultane ous outburst in the opt ical regime took 
place around 1997 (jTornikoski et al 19991) . Even in the 7- 
ray regime the source has been detected by the telescopes 
EGRET on-board CGRO and Fermi with a luminosity , 
L 7 = 5 x 10 47 erg/s dNolan et al.lll993l : lAbdo et alJl200flli . 

Within the framework of the mojaveQ program 
(Monitoring of Jets in Active galactic nuclei with VLBA 
Exp eriments) CTA 10 2 has been monitored since mid 1995 
(see lLister et al.1120091 and Lister et al. (these proceedings, 
p. 159)). The results of these intensive observations deliver 
a detailed picture of the morphology and kinematics of 
this source. The results of the kinematic analysis show an 
apparent velocity of the features in the jet between 0.7 c 



and 15.40 c ((Lister et al.ll2009h . A multi-frequency VLBI 
study inclu ding data at 90 GHz, 43 GHz and 22 GHz was 
reported bv lRantakvro et al.l (|2003l) . The results from the 
multi-frequency VLBI observations were combined with 
the continuum monitoring performed at single dish obser- 
vatories at 22 GHz, 37 GHz, 90 GHz, and 230 GHz. Within 
this multi-frequency data set (November 1992 until June 
1998) CTA 102 showed a major flare around 1997. The 
authors could conclude that this flare was connected to 
the ejection of a new jet feature. The inferred apparent 
speed of 11c, combined with the frequent and rapid flar- 
ing events throughout the electromagnetic spectrum, lead 
to the picture of a highly relati vistic jet. Thi s pictu re is 
supported by t he con clusions of Jorstad et al. (2005) and 
iHovatta et al. ( 20091 ). They found Lorentz factors, L, of 
17 and 15, and respectively, Doppler factors, S, between 
15 and 22. 

In April 2006, CTA 102 underwent a major radio flux 
outburst and we will present the result of our analysis 
of this flaring event using single dish and multi-frequency 
VLBI observations. We found indications that this flar- 
ing event was created by the interaction between a re- 
collimation shock and a traveling shock wave. 

2. Observations 

2.1. Single-Dish Light Curves 

In our analysis we concentrated on the radio flare around 
April 2006 and used single dish observations spanning 
from 4.8 GHz to 340 GHz (see Fig.[Tj). The flare is clearly 
visible at all the frequencies with increasing time delays 
towards smaller frequencies. The highest flux density of 
about 10 Jy is measured at 37 GHz. 

A self-absorbed synchrotron spectrum is described by 
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Fig. 1. Single dish light curves for CTA 102, centered 
around the 2006 radio flare. Used telescopes are 4.8- 
14.5 GHz UMRAO (M. F. Aller and H. D. Aller ), 
37.0 GHz Metsahovi (A. Lahteenmaki) and 230-340 GHz 
SMA (M. Gurwell) 
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Fig. 2. The 2006 radio flare in the turnover frequency 
- turnover flux density plane. The numbers indicate the 
temporal evolution of the flare. 



where S (v) is the flux density, v\ is the frequency at which 
the opacity t s = 1, and at and «o ar e the spectral in- 
dices for the optically thick and optically thin parts of 
the spectrum, respectively. The turnover frequency, v m , 
and the turnover flux density, S m , can be calculated from 
the first and the second derivative of synchrotron spec- 
trum and they can be regarded as the characteristics of 
the spectrum. 

For the spectral analysis we interpolated the data to 
a time sampling of At = 0.05 yr and subtracted an un- 
derlying quiescent spectrum. This homogenous spectrum 
(v m = 1.70 GHz, S m = 4.08 Jy, a = -0.35 and a t = 5/2) 
was created from archival data. The spectral evolution 
of the 2006 radio flare in CTA 102 is presented in the 
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Fig. 3. Evolution of the spectral parameters along the jet 
axis derived from the 2005.39 multi-frequency observa- 
tions of CTA 102. The estimate of the uncertainties are 
difficult due their dependence on the image alignment and 
the uneven uv-covcragc between the different frequencies. 
Note that the individual pixels are not independent. An 
estimate for the errors, assuming correct alignment, on the 
turnover frequency and flux density is around 15% of the 
derived value (see Discussion). 



turnover frequency - turnover flux density (v m — S m ) plane 
(see Fig.©. 

The flare starts around 2005.6 at a high turnover fre- 
quency (u m ~ 300 GHz) and low turnover flux density 
(Sm ~ 3 Jy). During the first 0.3 yr the turnover flux den- 
sity, 5 m , is increasing (to S m ~ 6.5 Jy) while the turnover 
frequency, v m is d e creasi ng (to v m ~ 120 GHz). Following 



Marscher fc Gear (1985) we could identify this stage as 



the Compton stage, where Compton losses are the domi- 
nant energy loss mechanism. The next stage in the shock- 
in-jet model should be the synchrotron one. This stage is 
characterized by a less prominent changes in the turnover 
flux density, S m , while the turnover frequency, v m is still 
decreasing. One could consider the time between 2005.8 
and 2005.9 as a possible candidate for this stage. During 
this time span of 0.1 yr the turnover flux density is slightly 
increasing (to S m ~ 7.0 Jy) while the turnover frequency 
keeps on decreasing (to u m ~ 90 GHz) . In the final stage 
the energy losses are dominated by the adiabatic expan- 
sion of the jet and the relativistic shock, the adiabatic 
loss stage. During this stage the turnover flux density and 
and the turnover frequency are decreasing. The adiabatic 
losses start to dominate the spectral evolution between 
2005.9 and 2006.0. 

The increase of the turnover flux density starting from 
2006.0 on and reaching a peak value of S m ~ 8.5 Jy in 
2006.3 can not be explained by the shock-in-jet model. 
After 2006.3 the turnover flux density decreases with de- 
creasing turnover frequency. A detailed analysis of the sin- 
gle dish observations will be presented elsewhere. 
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Fig. 4. 15 GHz VLBI image of CTA 102 observed on 6th 
of January 2007. 



2.2. Multi-Frequency VLBI Observations 

Since single-dish observations do not provide structural 
information of the jet we used as well multi-frequency 
VLBI observations covering a frequency range from 2 to 
86 GHz for accessing the neighborhood of the AGN central 
engine. CTA 102 was observed with the National Radio 
Astronomy Observatory's Very Long Baseline Array at 
three epochs (March 9th 2005, April 14th 2006 and 8th 
June 2006) using the all 10 antennas of the array. After 
calibration of the raw data, using the standard AIPS pro- 
cedures and model fitting in DIFMAP, we performed a 
core-shift analysis (Ar oc v 1 / kr ) following lLobanovl (1998) 
and a spectral analysis on the pixels along the jet axis 
using Equation [T] From the 2005.39 multi-frequency ob- 
servations of CTA 102 we concluded from the core-shift 
results that the source is in equipartition (jet particle en- 
ergy density equals magnetic field energy density). This 
circumstance is reflected by a value of k r — 0.98 ± 0.03 
and leads to an absolute distance from the central ob- 
ject Tcore^gGHz = 3.29 ± 1.12 pc and to a magnetic field 
B a 

ues (especially is B 
B oc vtSZ 2 . 



>GHz = 0.40±0.14 G. From the derived spectral val- 
S m ) we calculated the magnetic field 



Figure O presents the evolution of the physical param- 
eters deduced by our analysis as a function of the distance 
from the core. The results show an increase in the turnover 
frequency, u m , (upper right panel), the optically thin spec- 
tral index, ao, (lower right panel) and the magnetic field, 
B, at a distance of 0.2 mas away from the core. This be- 
havior could be an indication for an re-collimation shock 
(standing shock wave) at this position. 
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Fig. 5. Properties of the C12 component from the 
MOJAVE observations. Left: Separation from the core in- 
clusive linear fit to the acceleration event (dashed line). 
Right: Evolution of the flux density and time of beginning 
of the acceleration event (dashed line, see text for explana- 
tion). The uncertainties for the flux densities are around 
5% of the total intensities and the positional errors are 
of 1/5 of the restoring beam dimensions (typically around 
0.01 mas) 

2.3. MOJAVE Observations 

To check the evidence for a re-collimation shock at a dis- 
tance of 0.2 mas from the core we used the kinematic 
analysis of t he 15 GHz VLBI monitoring of CTA 102 



( Lister et al. 20091 ). A high resolution VLBI image of 



CTA 102 observed at 15 GHz on January 2007, showing 
the location and size all fitted circular Gaussian compo- 
nents, is presented in Fig. 2) The kinematics and the evo- 
lution of the flux density of a fitted component labeled as 
C12 could be interpreted as the interaction of a moving 
shock wave with the re-collimation shock: the separation 
from the core (r ~ 0.2 mas) remains constant until mid 
2006, when an acceleration event takes place together with 
a sharp drop in the flux density (see Fig. [5]). 

By fitting the acceleration part of the C12 trajectory 
and correcting for the core-shift we derived a value for 
the ejection time, t c - } = 2005.83 ± 0.05 yr. This time cor- 
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responds to the first peak in the v m -S m plane. The sec- 
ond peak is located slightly before the time of the be- 
ginning of the acceleration event (possible intersection 
of the moving shock wave with the re-collimation shock) 
frntersec ~ 2006.6 yr. Using a viewing angle of $ = 2. 6 ±0.5' 



(Jorstad et al 



2005) we derived a velocity /3 = v/c 
0.998 ± 0.042 and a bulk Lorentz factor T = 17.7 ± 0.7. 



3. Discussion 

The evolution of the 2006 radio flare in the Vm-S m plane 
could be explained by the interaction of a moving shock 
wave with a re-collimation shock. Re-collimation shocks 
are stationary features in non-pressure matched jets and 
lead to local increase in pressure and ch ange in the orien- 
tation and value of the magnetic field (|Dalv k, Marscher 
19881: iFallel 119911 IPerucho fe Marti 120071) . This behavior 
can be seen in the increase of the turnover frequency, v m 
and the increase in the magnetic field, B at a distance of 
0.2 mas in the evolution of the spectral parameters derived 
from multi- frequency VLBI observations (see Fig. [3J . The 
first peak in the v m -S m plane corresponds to the ejec- 
tion of the traveling shock wave around i c j ~ 2005.85 yr. 
The trajectory of this new feature can not be observed 
at 15 GHz due to the limited resolution. The interaction 
of the moving shock wave with re-collimation shock takes 
place around mid 2006 which could lead to the second 
peak in the t'm-S'm due to shock acceleration in a region 
of increased pressure and magnetic field. During the inter- 
action of the two waves the re-collimation shock is pulled 
away by the moving one and should appear after some 
time again at the same position. One possible explanation 
why this reaction is not detected in the 15 GHz kinematics 
could be a rarefaction wave. This wave is traveling in the 
wake of the moving shock wave and decreases the pressure 
behind the shock front. Together with the limited dynam- 
ical range of the VLBI this could be the reason for the 
non-detection of the re-collimation shock after the colli- 
sion with the traveling one. 



4. Conclusions and Outlook 

The combination of single-dish observations with multi- 
frequency, densely time-sampled VLBI monitoring is a 
powerful approach which can contribute towards a bet- 
ter understanding of flaring events. Using a high-quality 
set of observational data we presented a possible scenario 
which could explain the observed flare by a shock-shock 
interaction and derived estimates for the physical param- 
eters of the jet and the traveling shock wave. 

The spectral values presented have been derived from 
a spectral analysis applied on multi-frequency VLBI ob- 
servations. This technique is sensitive to the correct im- 
age alignment and to effects of the uneven uv-coverage 
between the frequencies. The effect on the derived param- 
eters due to small misalignments could lead to significant 
changes in the spectral values, especially in the magnetic 
field (B oc v^Sm)- Further analysis of the influence of the 



alignment is needed to provide adequate error bars for the 
derived values. 

Besides quantifying the uncertainties on the observa- 
tional parameters we started to test our assumption of 
the shock-shock interaction using 2D relativistic magneto- 
hydrodynamic simulations. These simulations will help 
us to understand the formation of re-collimation shocks 
in magnetized jets and their interactions with traveling 
shocks. 

It is also of interest to investigate possible correlations 
between 7-ray flares and the collision of re-collimation 
shocks and traveling shock waves. These investigations 
could help to clarify the question where in the jet the 
high energy radiation is generated. 
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